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Abstract. The fundamental notions of radiative transfer, e.g., Lam- 
bert's cosine rule, are studied from the point of view of flux and stress 
theory of continuum mechanics. For the classical case, where the ra- 
diance is distributed regularly over the unit sphere, it is shown that 
Lambert's rule follows from a balance law for the transfer of radiative 
power in each direction u of the sphere, together with the appropriate 
Cauchy postulates and the additional assumption that the correspond- 
ing flux vector field i„ be parallel to u. An analogous theory is presented 
for the irregular case where the distribution of radiance on the sphere is 
given as a Borel measure. 



1. Introduction 

We study the basic elements of radiometry from the point of view of flux 
and stress theory of continuum mechanics. Specifically, let x be a point on 
the boundary dR of some region R in space, let u be some unit vector and let 
6 denote the angle between u and the unit normal n to dR at x. Then, we 
recall that ignoring the dependence on the wavelength, the radiative energy 
flux density, the irradiance, E, out of dR at x is traditionally given in terms 
of the radiance field (radiative intensity) I(x, u) by 

E = [ I(x,u) cos 6{u)du(u), (1.1) 
Js 2 

where S 2 is the unit sphere containing all directions u and u denotes the 
solid angle on the sphere. The relation 

dE 

— = I(x, u) cos 9 (1.2) 

diO 

is usually referred to as Lambert's cosine rule {e.g., |Pre65[ p. 20] |SH92l p. 
20]). Thus, the total radiative energy flux out of R is given by 

P= / I(x,u)cosO(u)du)(u) 

JdR Us 2 

We show below that these traditional relations and more general expres- 
sions for the radiative energy flux follow from the theory of Cauchy fluxes 



dA (1.3) 
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supplemented with a single additional postulate. For the case where the radi- 
ance distribution at a point is a real valued function defined on S 2 , Cauchy's 
postulates are applied to the flux of energy in a generic direction u to yield 
a flux vector field \ u . Alternatively, and this is the approach used later for 
the more general situation where the distribution at a point may be as sin- 
gular as a measure, the collection of flux vector fields, {i u }, u £ S 2 , may 
be viewed as an infinite dimensional stress-like tensor field in the following 
sense. Given a point x and the unit normal vector n as above, the tradi- 
tional stress tensor determines the 3-dimensional vector of traction on dR 
at x. In analogy, the infinite dimensional tensor corresponding to radiance 
determines the distribution of flux on the sphere S 2 at x, an element of 
an infinite dimensional vector space C°(S 2 ) of continuous functions on the 
sphere. From a different point of view, the situation is analogous to mixture 
theory where a flux vector field is associated with each of the constituents. 
It is noted that the constituents of the analogous mixture should be labeled 
by a continuous variable on the sphere. Thus, the approach presented here 
offers a local point of view to the apparent non-local character of radiom- 
etry. One should also mention that a formal derivation of radiometry can 
be based on the principles of electromagnetism (e.g., [BW991 [Ger39]), where 
the radiance is obtained as average of the Poynting vector. 

The flux fields {i u } have an additional property that is characteristic of 
radiation: i u is parallel to u. This additional property implies the traditional 
law of radiometry as presented above. 

Replacing the space of continuous functions C°(S 2 ) by the space M(S 2 ) 
of Borel measures on the sphere, enables one to consider irregular radiation 
distributions such as monodirectional radiation. The property of radiation 
as described above cannot be applied directly to irregular distributions of 
radiation on the sphere because i u is not well defined. However, an analogous 
generalized formulation of the property is suggested by using some basic 
notions of measure theory. 



2.1. Traditional Cauchy Fluxes. Traditionally, Cauchy's flux theory con- 
siders an extensive scalar property p of regions R in the physical space which 
is represented here by R 3 . It is assumed that for each region R C R 3 , there 
is a field ipR : dR —> R, the flux density of p out of R, that represents the 
total flux of p through the boundary dR in the form 



A collection {i/jr} of fields on the boundaries of all regions R C R 3 will be 
referred to as a flux system. Cauchy's flux theory is developed on the basis 
of the following assumptions. 



2. Radiance Flux Fields: Classical Analysis 
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Boundedness (Balance). Let \R\ denote the volume of the region R. It is 
assumed that there is a positive constant C such that 

(2.2) 

The boundedness assumption above is usually motivated by a balance 
principle. It is assumed that the total flux of the property p out of R is 
equal to the rate of production of the property within R minus the rate of 
change of the total of p in R. Thus, let p be the density of p in IR 3 and let s 
be the density of the source of p in IR 3 , then 

f pdV + [ ipR dA = f sdV. (2.3) 

JR JdR JR 

The last balance equation, once supplemented by appropriate boundedness 
assumptions for p and s, will imply (I2.2p . 

It is observed that this assumption rules out irregular sources of the prop- 
erty such as those concentrated on surfaces, lines, or points. 

Cauchy's postulate of locality. Considering the dependence of the field ipR 
on the region R, Cauchy's locality postulate implies that this dependence is 
of a very short range. Specifically, it states that for x £ dR, iPr{x) depends 
on R only through the unit normal n(x) to dR at x. Let S 2 denote the 
2-dimensional unit sphere in IR 3 containing the collection of unit vectors. It 
follows that there is a function r : IR 3 x S 2 ->■ R such that 

ip R (x)=T(x,n(x)). (2.4) 

Regularity. It is assumed that r is a smooth function. 

Naturally, we will refer to a flux system satisfying these assumptions as a 
Cauchy Flux System. 

Cauchy 's Flux Theorem. The boundedness and regularity assumptions above 
imply that the dependence of r on n is linear. Since any linear function of n 
may be represented by an inner product with n, it follows that for a Cauchy 
flux system there is a vector field T : IR 3 — > IR 3 such that 

if) R {x) = t(x, n(x)) = T(x) ■ n(x). (2.5) 

We refer to T as the flux vector field associated with the property p. 

Consider for example the case where the property under consideration is 
the radiation energy so that ^ r is interpreted at the total flux of radiation 
energy through the boundary of a region R. Then, under the foregoing 
assumptions, Cauchy's flux theorem implies that there is a vector field q 
such that the flux of radiation energy out of the boundary of a region R is 
given by (12. 5|) . In fact, for this case, q ■ n is the irradiation as implied by 
equation (ll.3p . 
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The Differential Balance Equation. Using Cauchy's theorem in the balance 
equation (12.3|) we have 

pdV + [ T-ndA = [ sdV. (2.6) 
: JdR Jr. 

Using Gauss's theorem one concludes that 

div T + p = s in fl, 

(2 7) 

T(x) ■ n(x) = ip R (x) on dR. { ' J 

Remark 2.1. Let {i^r} be a flux system. Assume that there is a differentiable 
vector field T such that for each region R, tpR = T ■ n, then, 

y R = f ip R dA= [ T-ndA = [ divTdU. (2.8) 

JdR JdR JR 

It follows that if divT is bounded, the boundedness assumption (12. 2p holds. 
This observation is an inverse to Cauchy's theorem. 

For detailed, generalized and technical presentations of Cauchy's flux the- 
ory see for example |Sch07| and references cited therein. 

2.2. Radiance Systems. The space of continuous real valued mappings 
defined on the sphere will be denoted by U. For a regular region R in space, 
a radiance field over OR is a mapping 

Kr : dR — > U. (2.9) 

For x G dR and u € S 2 , Kr(x)(u), is interpreted as the density of radiation 
power per unit area at x € dR crossing dR in the direction u. It should be 
emphasized that this interpretation is given here without any justification 
in order that the presentation is adapted to the context of radiation (see 
Section E3J. 

Clearly, the radiance field may be regarded as a function K' R : dR x 
S 2 — > K by K' R (x,u) = Kr(x)(u). We prefer using Kr because of the 
generalization of Section [3l 

For a given u G S* 2 , the radiance field at the direction u is the function 

K RtU : R 3 — > R, K RjU (x) = K R (x)(u). (2.10) 

The total emitted power density at x, the irradiance or radiant emittance is 
therefore given by 



M%)= / K R (x)(u)doj(u) = / K R (x)(u)doj(u), (2.11) 
Js 2 Js 2 

where oj is the solid angle measure. The total power emitted from dR is 

Pr= f il> R (x)dA= [ ([ K R (x)(u)dLj(u) ) dA (2.12) 

J dR JdR \JS 2 / 

Fubini's theorem implies that we can write 

P R = f ( [ K R (x)(u)dA(x)]dcj= [ PR, u du(u), (2.13) 

JS 2 \JdR J JS 2 
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where 

Pr. u = [ K R (x)(u)dA(x) = ( K R , u {x)dA{x) (2.14) 

JdR JdR 

is the total power emitted from R in the direction u £ S 2 {directional emissive 
power m |Mah02| ). It is noted that for a fixed aeS 2 , P R u is the total flux 
of the density K R>u (x) = K R (x)(u) — a scalar flux density over dR. 

A radiance system is a collection {K R } for all regions R C R 3 . The basic 
objective of flux theory is to study the dependence of the radiance K R on 
the region R. 

We will say that a radiance system {K R } is a Cauchy radiance system if 
for each u £ S 2 , the scalar valued flux system {K R ^ U } satisfies Cauchy's 
postulates. From Cauchy's theorem for scalar fluxes, it follows that for 
a Cauchy radiance system, for every u £ S 2 , there is a vector field, the 
radiance vector field for the direction u, 



! , such that i u (x) ■ n(x) = K R ^ u (x), 



(2.15) 



where n{x) is the unit normal to dR at x. 

Using the radiance vector field, Equation (I2.14P may be written as 



JdR 



i„ ■ n dA. 



(2.16) 



The existence of the radiance vector field for the direction u enables one 
to define the field i over I 3 x I 3 , whose value is a function over the sphere 

by 

i(n, x)(u) = i u (x) ■ n. (2-17) 
3 — > U as a radiation density tensor and the 



We will refer to i : R x 
expression for the power becomes 



Pr 



or 



r)R 



S 2 



{i u (x) ■ n(x))duj(u) 



i(n(x), x)(u) dco(u) 



dA, 
dA, 



(2.18) 



Clearly, i is linear in n and using e,- to denote the j'-th standard base 
vector in R 3 , one may define 



so that 



ij(x) = i(ej,x), 
i(n,x) = ^2ij(x)rij, 



(2.19) 
(2.20) 
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and 



Pr = f dR \J2 i i( x ) n i) W duj ( u ) 
= / E , W LA ' 

JdR, 



dA, 



d.4, 



(2.21) 



where 



<?j = / ih{x)){u)duj(u) 
Js 2 

are the components of the total radiation energy flux vector field q. 



(2.22) 



Remark 2.2. It is noted that the equations above make it possible to obtain 
the total flux across surfaces that are not necessarily the boundaries of re- 
gions. For any smooth surface S, one has to replace the integration over dR 
by integration of S. 

2.3. The Source Term and Differential Balance Equation. Using Gauss's 
theorem in Equation (12.16P one has 



Pr,u = I div i u dV 
Ir 



(2.23) 



and 



Pr 



div i„ dV 



R 



div \ u du 



du, 
dV. 



(2.24) 



A balance equation for the energy flux in the direction u G S 2 will be of 
the form 



I p u dV+ [ K R>u dA= [ 

JR JdR JE 



S„. dV 



(2.25) 



R 



and with Equation (I2.15|) 

p u dV + I \ u -ndA 



s„.dV. 



IR JdR JR 

The corresponding differential balance equation is 

divlu + p u = s u . 

Integration over the unit sphere will give the balance equation for the total 
energy. 



(2.26) 
(2.27) 
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2.4. The Basic Assumption for Radiance. It is noted that thus far very 
little properties of radiation were used. One could replace the sphere S 2 by 
some other space B so that the space U would be the collection of continuous 
functions on B. For example, if the space B is the set {1,2,3}, the space 
U of real valued functions on it is identical to R 3 by identifying Wj with 
w(i) for each i € B = {1,2,3}. Clearly, continuity is insignificant in this 
example. Such a construction will lead us to traditional stress theory as we 
remark below. In fact, if we do not insist that the functions considered be 
real valued, the set U can be any set. We could consider a system {K^ u }, 
for each region R and each u £ B and require that it satisfies Cauchy's 
postulates for and fixed u. For each u £ B, Cauchy's theorem would imply 
the existence of a flux vector field \ u . 

Among such theories, what characterizes radiance theory is the assump- 
tion that radiance vector field i n is parallel to u, i.e., 



where i u = \i u \ is a positive real valued function on R 3 . 
Equation (|2.15p that 

K(x, n, u) = i u (x)u ■ n. 
Denoting by 6 the angle between the vectors u and n, 

K(x, n, u) = i u (x) cos 6. 



(2.28) 
It follows from 

(2.29) 
(2.30) 



With this basic assumption, the interpretation of Kr(x)(u) as the flux of 
energy flowing in the direction u is justified. For vectors u that point out 
of R, Kr(x)(u) is the power flux density outgoing from R and for vectors 
u that point into R, Kr(x)(u) is the incoming power density. One can 
also consider surfaces that need not be the boundaries of regions so that 
"incoming" and "outgoing" are meaningless. In such circumstances Kr(x)(u) 
and Kr(x)(—u) are specified independently in order to distinguish between 
fields such as those generated in the middle between two equal light sources 
and fields where there are no light sources. 

Equation (12. 16ft may now be written as 



Pr,' 



i 7 ,u ■ n dA. 



(2.31) 



Thus, the expression (I2.18P for the power becomes 



Pr 



dR 



dR 



S- 



i u (x) cos 6(u) du(u) 



dA, 



i u (x)u duj(u) 



(2.32) 



n(x) dA, 



so that the radiation energy flux vector field is given by 



q{x) 



i u (x)u duj(u). 



(2.33) 
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Alternatively, we may write 



Pr 



OR 

S u ■ u du, 



t (x)n dA 



uduj(u), 



where 



r)R 



t (x)ndA 



(2.34) 



(2.35) 



represents the total contribution of the radiation in the direction of u (cf. 
radiant intensity J? in [AK961 Equation (2.1.9)]. 

Comparing Equation (|2.32p with traditional expositions of radiation the- 
ory e.g., |Pla591 ICha60l IAK961 IMod93l ISob63| . it is noted that i u (x) is the 
scalar radiance density e.g., K(x,u), that they use. 

Remark 2.3. In order to get some insight to the various objects defined and 
derived here, we compare the notions under consideration to stress theory of 
continuum mechanics. In continuum mechanics, the traction vector t, deter- 
mined by its 3 components ij, depends on the point x and external unit nor- 
mal n. Thus, while t is determined by the 3 components ij the radiation field 
is determined but the infinite number of "components" Kr <u (x) = Kr(x)(u) 
parametrized by the direction u. 

The vector field i u that induces the radiance, is thus analogous to the i-th 
row of the stress tensor a, Ui, that gives the i-th component of the traction 
vector by 

ti = Oi ■ n. (2.36) 

Again, just as there are 3 such rows in the stress tensor, there are uncountably 
infinite row vectors i u , each for any value of u. 

The property of radiation implying that the radiance flux vector for the 
direction u is parallel to u is therefore analogous to the property that the 
i-th row of the stress tensor is parallel to the i-th coordinate axis. Thus, the 
analogous property for stresses is that the stress matrix is diagonal. 



We may now apply Green's theorem to Equation (|2,35p to obtain 



Vz„ dV. 



(2.37) 



R 



Thus, the total flux out of R may be expressed as 



Pr 



Vi„ dV 



■ u du, 



Vi 7/ • u duo 



R US 2 



(2.38) 



dV. 



This also follows from 



divi u = div {i u u) = Vi u • u. 
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Evidently, the same result could be obtained by applying Green's theorem 
to the second line of Equation (12. 32ft . The differential balance equation for 
the radiation in the direction u, Equation (|2.27p . assumes the form 

Vi u ■ u + p u = s u . (2.39) 

Remark 2.4. Assume that p u = s u = in the last equation. If follows that 

Vi u ■ u = 0, (2.40) 

i.e., the directional derivative in the direction u of the radiance in the di- 
rection u vanishes. This result, formulated as the conservation of i u under 
translation in the direction u, is sometimes referred to as the fundamental 
theorem of radiometry (e.g., |Boy83 pp. 18-19]). 



For the differential balance equation (12. 7ft corresponding to the total ra- 
diation energy 

divq + p = s, (2.41) 

one has 

divq = / Vi u -uduj. (2.42) 
Js 2 



Remark 2.5. From Remark 1 2. 1| it follows that if one starts from the standard 
expression (|2.32p . then, Cauchy postulates follow for different iable fields i u . 
In other words, the Cauchy postulates are both necessary and sufficient for 
the traditional radiation expression. Not merely some, but all differentiable 
radiance fields represent Cauchy radiance systems. Specifically, it follows 
from Equation (|2.3ip that 

Pr, u = [ div(i u u)dV (2.43) 
Jr 

and so 

|Pr,„| < max|div(i u w)(x)| \R\ , (2.44) 



the boundedness expression for Pr 



2.5. Virtual Power. Let u be a fixed direction is space and let w u : dR — > 
R be a field. One may consider the action 

Pr,u( w u)= K R<u (x)w u (x)dA(x). (2.45) 

JdR 

We interpret Pr )U (w u ) as follows. Let W u : dR — s> M + be a differentiable 
function that we interpret as a distribution of radiance meters over the 
boundary that measure the radiance in the direction u so that W u (x) is 
the density of the meters at x. The field w u is conceived as a variation of 
W u . Thus, while W u (x) should be positive by our interpretation, this limita- 
tion does not apply to w u (x). Hence, Pr )U (w u ), is interpreted as the change 
of the total power of radiation propagating in the direction u measured by 
the distribution of meters under the variation w u of the distribution. For 
short, we will refer to Pr, u {w u ) as the virtual power. 



10 



REUVEN SEGEV AND JOE GODDARD 



Using the radiance vector field i u we may write the virtual power as 

Pr,u(w u ) = / w u (x)i u (x) ■ n(x)dA(x). (2.46) 
J OR 

The last equation may be transformed using Gauss's theorem into 

Pr, u (w u ) = / div(to„iu) clV, 
JR 



Vw u ■ i u dV + / w u div i u dV, 
R JR 



(2.47) 



and with Equation (|2,27p we obtain 

/ K R)U w u dA+ / w u {p u -s u )dV= / Vw u -i u dV. (2.48) 

J OR JR JR 

Using the basic assumption for radiance (|2.28p we finally arrive at 

/ K R)U w u dA+ / w u (p u -s u )dV= / i u Vw u ■ udV, (2.49) 
JdR Jr jr 

where evidently Vw u • u is the directional derivative of w u in the direction 



of u. 



3. Measure Valued Radiance 



3.1. Preliminaries. In the foregoing analysis, the flux density of radiation 
energy is distributed continuously on the sphere S 2 . For example, for a 
point x G dR, a function K x = Kr{x) : S 2 — > R describes the distribution of 
radiation power, where for u € S 2 , K x (u) was interpreted as the radiation 
density emitted in the direction of u. It is noted that if one wishes to 
consider power radiated in one particular direction Uq, a distribution such 
as K x above is not general enough. One has to include mathematical objects 
such as the Dirac delta distribution, 5 Uo , representing the emission of power 
in the particular direction uq and which is a measure on the sphere. It 
is therefore assumed henceforth that the distribution of radiation, or the 
radiance, at x € dR, is a Borel measure on the unit sphere. Thus, if we 
denote the radiance at x by J x , for any Borel subset D of the sphere S 2 , 
J X {D) denotes the measure of D indicating the total power density per unit 
area emitted through the pencil at x subtended by D. In particular, the 
irradiance is given by 



E = J X (S 2 ) = I dJ x . (3.1) 
Js 2 



is 2 

In the traditional continuous case, the measure J x is absolutely continuous 
relative to the solid angle measure, or area measure on the unit sphere, so 
that there is a Radon-Nikodym derivative 

K x = ^, K X :S 2 ^R (3.2) 
dui 
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and 

J X (D)= I K x du. (3.3) 



D 

(See for example |Pre65| IArv95| who postulate the continuity requirement 
right from the start.) As before, the function K x {u) indicates the radiance 
at x out of R in the direction of u and the total density of power, emitted 
at x in all directions is 

/ K x du. (3.4) 
Js 2 

We will denote by M the vector space of Borel measures on the unit sphere, 
and so, J x £ M (in comparison with K x € U). The vector space operations 
are naturally defined in M by (aJ x + a' J' X ){D) = aJ x (D) + a'J' x (D). In 
addition, one can define a norm on M providing it with the structure of a 
Banach space. 

It is recalled that the space of Borel measures M may be identified with 
the dual space U* of the space U of continuous real valued mappings on the 
sphere. Specifically, the measure J x may be regarded as a continuous and 
linear functional J' x : U — >• M. such that " 



J' x (v)= / v(u)dJ x (u). (3.5) 
Js 2 

The expression above is interpreted as the virtual power density at x for the 
change of distribution of radiance meters v at x. In the sequel, we will not 
distinguish in the notation between the measure and the linear functional it 
induces. In terms of this action, a natural norm on M is given by 

II 7 II J *( V ) /Q fiN 

\\J X \\ = sup n ,. , (3.6) 

veu IMI 

where = max ugS 2 

3.2. Measure Valued Radiance Systems and Radiance Tensors. In 

analogy with the definition of the radiance field over dR in (|2.9p , the basic 
assumption is that for every region Eci 3 there is a smooth function, the 
measure valued radiance, 

J R :dR — ► M. (3.7) 
The collection {Jr}, for all regions R is a measure valued radiance system. 

For a given measure valued radiance system and an element v € U, let 
{Jr,v} be the flux system such that 



Jr,v(x) = Jr{x){v) = / vd(J R (x)). (3.8) 
Js 2 

i.e., the action of Jr(x) on v is the action of a continuous linear functional 
on a continuous function defined on the sphere. Similarly, one can evaluate 
the integral of the measure Jr(x) over any integrable Borel subset D C S 2 
to yield 

J R (x)(D) = f d(J R (x)). (3.9) 

JD 
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In analogy with Section T2.2I we now assume that for each fixed v G U, the 
measure valued radiance systems {Jr, v } satisfies Cauchy's postulates. Let 
v G U be a fixed distribution. It follows from Cauchy's theorem for scalar 
valued flux systems that there is a vector field I„ = £^ Ivi e i : K 3 — )• R 3 such 
that 

Jr,v(x) = Jr(x)(v) = I v (x) ■ n{x) = ^l v i(x)rii(x). (3.10) 

i 

for any x G dR. It follows from Equation (|3.10p that I v and I v i depend on 
v G [/ continuously and linearly. As I«(x) are real, there are three elements 
Ii(x) G M such that I v i(x) = li(x)(v). Hence, 



Jr,v(x) 



J2ni(x)Ux){v) =J2Mx) I vd(Ux)). (3.11 



Clearly, G M are the three flux distributions at x corresponding to the 

area elements which are perpendicular to the standard base vectors in M 3 . 

Using the notation L(R 3 ,M) for the space of linear mappings M 3 — > M 
(which are also continuous as they are defined on a finite dimensional space) , 
one can now define the mapping 

I:M 3 — ► L(IR 3 , M) = M (g) R 3 (3.12) 

by 

(I(x)(n))(v) = I v (x) n = ^2 ni(x)Ii(x)(v). (3.13) 

i 

Observe that for any net 3 , 

^2li(x)rii = I(x)(n), 

i 

= Ux)(^2nieA, (3.14) 

^ i ' 

= ^2ml(x)(ei), 

i 

so that li(x) = I(x)(ej) G M. Thus, we may introduce the notation 

I = ^I i ®c <J (3.15) 

i 

and write 

J R = I - n = (^2*1 ® • (njej) = ^IjUj = I(n). (3.16) 
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The infinite dimensional tensor I will be referred to as the measure valued 
radiance tensor. In particular, for any Borel integrable set D, 



Jr(x)(D) = h(x)(D) ni (x) = I(x)(D) ■ n(x), 



i 



ni (x) / d(ii(x)). 



(3.17) 



D 



In the case where the measures I(x) may be represented by densities rel- 
ative to the solid angle measure cj on the sphere, one has 

dl(x) = i(n,x)doj (3.18) 

where i(n, x) is the radiance density tensor. 

Using the measure valued radiance tensor, the total energy flux out of a 
region R is given as 



Pr 



dR 



n(x) ■ d(I(x)){u) 



dA(x), 



dR 



si 



d(I(x))(u) 



n(x) dA(x), 



(3.19) 



dR 



I(x)(S 2 ) -n(x)dA(x), 



where, 



and Ii(x)(S 2 



I(x)(S 2 ) = y £l i (x)(S 2 )®e l 



(3.20) 



1,2,3, represent the total measures of the sphere. The 



energy flux vector field q is therefore given by 

q(x)= [ d(I(x))(u)=I(x)(S 2 ). 
Js 2 



(3.21) 



3.3. Measure Valued Radiance Tensors: Totals Approach. We present 
below an alternative derivation of the measure valued radiance tensor in 
terms of a total radiation distribution corresponding to a region R. In order 
to formulate the theory, one needs to use integration and differentiation of 
functions on ]R 3 valued in a Banach space — the space M in our case. (See 
for example |Lan69[ Chap. V] or |Die60[ Chap. VIII] for some details on 
calculus of functions valued in Banach spaces.) 

3.3.1. The Analogs of Cauchy's Postulates. As in the definition of the ra- 
diance field over dR in (13. 7p . the basic assumption is that for every re- 
gion R C M 3 there is a smooth function, the measure valued radiance, 
Jr : dR — > M. The total of the radiation corresponding to the region 
R is defined as the distribution 



Jr dA G M, 



dR 



(3.22) 
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where we integrate the measures Jr{x) over dR. In other words, for each 
measurable subset D C S 2 , &r(D) is the total flux of energy flowing in the 
directions included in the pencils determined by D from all points in dR. 

In case &r is absolutely continuous relative to the solid angle measure 
on the sphere, the directional emissive power Pr u is given by the Radon- 
Nikodym derivative 

d$p 

Pr, u = (3-23) 



The following assumptions are made in analogy with Section 12.11 

Locality: It is assumed that for a point x £ dR, Jr{x) depends on R through 
the outwards pointing normal n{x) to dR at x. 
It follows that there is a function 

J : R 3 x S 2 — YM, (3.24) 

to which we will refer as the Cauchy mapping, such that 

J{x,n{x)) = J R {x) (3.25) 

for any region R such that x 6 dR and n(x) is the outwards pointing normal 
to R at x. 

Using the Cauchy mapping we may rewrite the total as 

$ R = / J(x,n(x))dA. (3.26) 
JdR 

Boundedness (Balance): We assume that there is a positive number C such 
that 

<C|i2|, (3.27) 
where is the norm in M of the total. 

Regularity: It is assumed that the Cauchy mapping J is smooth. 

3.3.2. The Analog of Cauchy's Theorem. We outline a sketch of the proof of 
Cauchy's theorem for the current settings. Consider an infinitesimal tetra- 
hedron containing the point x £ ffi 3 such that for a = 0, 1,2,3, A a are the 
areas of the faces of the tetrahedron and n a are the unit normals to the 
faces. Thus, Equation (|3,26p assumes the form 

®R = Y1 Jfan«)Ax, (3.28) 

a 

and the boundedness assumption (|3.27p implies that 



\\®r\\ = 

Since n a A a = 0, we have 



J (x, Tloi)A a 



^C\R\. (3.29) 



n » = "I!^ n !" (3.30) 
P =i 
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which implies 




(3.31) 



As the size of the tetrahedron approaches zero, the right-hand side of the last 
equation tends to zero. Taking the limit, it follows that the left-hand side of 
the last equation vanishes. Since for any norm on M, \\Jo\\ = implies that 
Jo = € M, we obtain 

J U-ET n A =-E^ J (^)- ( 3 - 32 ) 
\ p=i / p=i 

The last equation implies that the dependence of J on its second argument 
n is linear. 

We conclude that there is a field 

I : M 3 — ► L(R 3 ,M) given by I(x)(n) = J(x 1 n). (3.33) 

Using the mapping I, one can write the analog of Cauchy's formula as 

I(x)(n) = Jr{x) (3.34) 

for every region R such that n is the outwards pointing normal to dR at x. 
Hence, I is identical to the measure valued radiance tensor of fj3. 13[) . 

Using the measure valued radiance tensor, the irradiance, the total of the 
radiation may be written as 

$r= I(x) • n(x) dA e M. (3.35) 
JdR 

Using the Gauss theorem for the measure valued field I we may rewrite 
Equation (13.351) as 

$ R = / divIdU. (3.36) 
Jr 

The measure may be integrated over the sphere (or any other Borel 
measurable subset thereof) to give the total power 

P R = f d$ R = [ d( f I(n)(x)dA (3.37) 

3.4. The Basic Assumption of Radiation for Measure- Valued Radi- 
ation Tensors. The basic assumption for radiance theory should be gener- 
alized so that it applies to the measure valued radiance tensor. The assump- 
tion cannot be formulated as in (I2.28P because in the general case one cannot 
assume that the vector density i u exists. Nevertheless, the basic constitutive 
assumption may be generalized as follows. 
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We recall that a measure fx, viewed as a continuous, linear functional 
acting on continuous functions, may be multiplied by an integrable function 
(f) to yield the measure <p fx defined by 

(<f>®H){f)=fi(<f>f), or equivalently, J f d(<f> y) = j f<t>d[i. (3.38) 

In addition, the Radon-Nykodim theorem implies that if there is a positive 
measure A such that \(D) = implies that fi(D) = 0, then there is an 
integrable function, the Radon-Nykodim derivative 

f-% (3-39) 

such that /x = / A. 

We now use a procedure as in |Edw65[ pp. 236-239]. For each x € ffi 3 
which is kept fixed in this paragraph, consider the positive measure 



|I(x)| = |I| (x) := 



, where, [!*(*)] (D) = pL,(x)(JD)] 2 . 

3 

(3.40) 

Clearly, for each j = 1,2,3, |Ij(x)(D)| |I(x) | for every measurable 
subset D. It follows from the Radon-Nikodym theorem that for each j there 
is an integrable function lj(x) defined on S 2 such that 

Ij{x) = Ij{x) |I(x)| and = 1 (3.41) 

j 

for almost all u G S 2 relative to the measure |I(x)|. For example, for any 
measurable subset D, 

[ d(I j (x))= [ Ij(x)d\l(x)\. (3.42) 
Jd Jd 

One may define the vector function I(x) : S 2 — >• IR 3 on the sphere by 

I(x)(u) = j;i i (x)(u) Cj (3.43) 

j 

and so 

I(x) = X)li(x) Ci = I(x) |I(x)| , (3.44) 

3 

where in the equation above indicates the product of the measure |I(x)| 
by the vector valued integrable function l(x). In other words, all the sin- 
gularities of the measures Ij(x) are included in |I(x)| and the "unit vector" 
I(x) distributes them to the various directions. 

The basic assumption of radiation theory may be formulated simply as 

I(x)(ti)=t*. (3.45) 
In other words, for each x € M 3 there is a scalar measure |I(x)| such that 

I(x) = lQ \1{x)\ , (3.46) 
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where 

l : S 2 -> S 2 (3.47) 

is the identity mapping on the sphere. Again, on the right-hand side of 
Equation (|3,46p is the product of a vector (or a sphere)-valued function and 
a measure. 

It follows from Equations (|3.13H3~T5|) that 

J(x,n) = I(x) ■ n = \I(x)\ {u ■ n). (3.48) 
Using Equation (|3.48p the total distribution may be written now as 

$ R = / |i| ( u ■ n) dA, (3.49) 

JdR 

and using Gauss's theorem 

®R = I y,lT^ @u i dV = I (V|I|)OtidV, (3.50) 
Jr j dx j Jr 

where the gradient of the measure valued function |I| is used as well as the 
notation (V |I|) u = £\ |^ Uj . 
For any measurable D C S 2 , 



[ d(I(x)-n) = [ d(|I(z)|«-n), 
Jd Jd 

= / u ■ nd |I(x)| . 
Jd 



(3.51) 



In the classical case where d |I(x)| = i u {x)duj we revert to 

/ d(l(x) ■ n) = / u ■ ni u {x)duj. (3.52) 
Jd Jd 
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